Studying Venom Toxin Variation Using Accurate Masses from Liquid Chromatography-Mass Spectrometry Coupled with Bioinformatic Tools.

This study provides a new methodology for the rapid analysis of numerous venom samples in an automated fashion. Here, we use LC-MS (Liquid Chromatography-Mass Spectrometry) for venom separation and toxin analysis at the accurate mass level combined with new in-house written bioinformatic scripts to obtain high-throughput results. This analytical methodology was validated using 31 venoms from all members of a monophyletic clade of Australian elapids: brown snakes (Pseudonaja spp.) and taipans (Oxyuranus spp.). In a previous study, we revealed extensive venom variation within this clade, but the data was manually processed and MS peaks were integrated into a time-consuming and labour-intensive approach. By comparing the manual approach to our new automated approach, we now present a faster and more efficient pipeline for analysing venom variation. Pooled venom separations with post-column toxin fractionations were performed for subsequent high-throughput venomics to obtain toxin IDs correlating to accurate masses for all fractionated toxins. This workflow adds another dimension to the field of venom analysis by providing opportunities to rapidly perform in-depth studies on venom variation. Our pipeline opens new possibilities for studying animal venoms as evolutionary model systems and investigating venom variation to aid in the development of better antivenoms.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes.

BACKGROUND: Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris). RESULTS: Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis. CONCLUSIONS: Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Prevalent bee venom genes evolved before the aculeate stinger and eusociality.

BACKGROUND: Venoms, which have evolved numerous times in animals, are ideal models of convergent trait evolution. However, detailed genomic studies of toxin-encoding genes exist for only a few animal groups. The hyper-diverse hymenopteran insects are the most speciose venomous clade, but investigation of the origin of their venom genes has been largely neglected. RESULTS: Utilizing a combination of genomic and proteo-transcriptomic data, we investigated the origin of 11 toxin genes in 29 published and 3 new hymenopteran genomes and compiled an up-to-date list of prevalent bee venom proteins. Observed patterns indicate that bee venom genes predominantly originate through single gene co-option with gene duplication contributing to subsequent diversification. CONCLUSIONS: Most Hymenoptera venom genes are shared by all members of the clade and only melittin and the new venom protein family anthophilin1 appear unique to the bee lineage. Most venom proteins thus predate the mega-radiation of hymenopterans and the evolution of the aculeate stinger.

Domain loss enabled evolution of novel functions in the snake three-finger toxin gene superfamily.

Three-finger toxins (3FTXs) are a functionally diverse family of toxins, apparently unique to venoms of caenophidian snakes. Although the ancestral function of 3FTXs is antagonism of nicotinic acetylcholine receptors, redundancy conferred by the accumulation of duplicate genes has facilitated extensive neofunctionalization, such that derived members of the family interact with a range of targets. 3FTXs are members of the LY6/UPAR family, but their non-toxin ancestor remains unknown. Combining traditional phylogenetic approaches, manual synteny analysis, and machine learning techniques (including AlphaFold2 and ProtT5), we have reconstructed a detailed evolutionary history of 3FTXs. We identify their immediate ancestor as a non-secretory LY6, unique to squamate reptiles, and propose that changes in molecular ecology resulting from loss of a membrane-anchoring domain and changes in gene expression, paved the way for the evolution of one of the most important families of snake toxins.

Diversity Begets Diversity When Diet Drives Snake Venom Evolution, but Evenness Rather Than Richness Is What Counts.

Snake venoms are primarily used to subjugate prey, and consequently, their evolution has been shown to be predominantly driven by diet-related selection pressure. Venoms tend to be more lethal to prey than non-prey species (except in cases of toxin resistance), prey-specific toxins have been identified, and preliminary work has demonstrated an association between the diversity of diet classes and that of toxicological activities of whole venom. However, venoms are complex mixtures of many toxins, and it remains unclear how toxin diversity is driven by diet. Prey-specific toxins do not encompass the molecular diversity of venoms, and whole venom effects could be driven by one, few, or all components, so the link between diet and venom diversity remains minimally understood. Here, we collated a database of venom composition and diet records and used a combination of phylogenetic comparative methods and two quantitative diversity indices to investigate whether and how diet diversity relates to the toxin diversity of snake venoms. We reveal that venom diversity is negatively related to diet diversity using Shannon's index but positively related using Simpson's index. Since Shannon's index predominantly considers the number of prey/toxins, whereas Simpson's index more strongly reflects evenness, we provide insights into how the diet-venom diversity link is driven. Specifically, species with low diet diversity tend to have venoms dominated by a few abundant (possibly specialised) toxin families, whereas species with diverse diets tend to 'hedge their bets' by having venoms with a more even composition of different toxin classes.

Assessing knowledge and awareness regarding snakebite and management of snakebite envenoming in healthcare workers and the general population: A systematic review and meta-analysis.

BACKGROUND: Snakebite envenoming is a serious and life-threatening medical condition that predominantly affects people living in rural communities across Africa, Asia, and Latin America. As our climate changes, there is a growing concern that negative human-snake interactions will increase. Our ability to prevent and manage snakebite requires effective antivenoms as well as knowledge regarding the prevention and management of snakebite among healthcare workers and affected communities across the globe. This systematic review aims to assess existing levels of knowledge regarding snakebite prevention and management in both healthcare workers and affected communities. METHODS: This review was conducted on studies reporting quantitative measurements to evaluate knowledge and practice regarding snakebite prevention and management published in major databases between 1 January 2000 and 31 December 2021. Random effects modelling was used to obtain the pooled proportion. Heterogeneity (I2) was tested, and sensitivity analyses performed. RESULTS: Out of 3,697 records, 16 studies from 12 countries assessing 7,640 participants were included. Four of the studies were ranked as good quality studies, 9 as fair, and 3 as poor. This study results demonstrated that 56% of the study population answered the knowledge question correctly (95% CI 48% to 63%, p < 0.001). High heterogeneity was observed (I2 = 97.29%), with marginal publication bias (Egger's regression test, p = 0.0814). Participants had relatively higher knowledge concerning use of antivenom as preferred treatment, followed by snakebite prevention, knowledge of signs and symptoms of snakebite, knowledge of first-aid, and knowledge of treatment. Participants had lower knowledge relating to types of snakes and the identification of snakes. CONCLUSION: Adequate knowledge about snakebites and its management among the general population and healthcare workers was 56%. Healthcare workers and communities across Asia showed higher relative knowledge compared to those in Africa and the Middle East. These data suggest that further education is needed in both the general population and among healthcare workers to ensure that appropriate preventative and patient management techniques are being utilised in snakebite endemic regions. Greater local awareness of the risks and appropriate management of snakebite is required to reduce the burden of snakebite mortality and morbidity.

Highly Evolvable: Investigating Interspecific and Intraspecific Venom Variation in Taipans (Oxyuranus spp.) and Brown Snakes (Pseudonaja spp.).

Snake venoms are complex mixtures of toxins that differ on interspecific (between species) and intraspecific (within species) levels. Whether venom variation within a group of closely related species is explained by the presence, absence and/or relative abundances of venom toxins remains largely unknown. Taipans (Oxyuranus spp.) and brown snakes (Pseudonaja spp.) represent medically relevant species of snakes across the Australasian region and provide an excellent model clade for studying interspecific and intraspecific venom variation. Using liquid chromatography with ultraviolet and mass spectrometry detection, we analyzed a total of 31 venoms covering all species of this monophyletic clade, including widespread localities. Our results reveal major interspecific and intraspecific venom variation in Oxyuranus and Pseudonaja species, partially corresponding with their geographical regions and phylogenetic relationships. This extensive venom variability is generated by a combination of the absence/presence and differential abundance of venom toxins. Our study highlights that venom systems can be highly dynamical on the interspecific and intraspecific levels and underscores that the rapid toxin evolvability potentially causes major impacts on neglected tropical snakebites.

Molecular mechanisms underlying the N-methyl-d-aspartate receptor antagonists: Highlighting their potential for transdiagnostic therapeutics.

The so-called "psychedelic renaissance" has stimulated expanded interest in several classes of drugs that appear to possess transdiagnostic effects in the treatment of mental health disorders, specifically. N-methyl-d-aspartate receptor (NMDAR) antagonists are one such class with diverse therapeutic potential. NMDARs mediate excitatory postsynaptic signalling in the central nervous system (CNS) and are integral to normal neurobiological processes including neuronal development, synaptic transmission, and plasticity, and thus involved in learning and memory. However, NMDAR hyper-function is also implicated in acute CNS trauma, neuropsychiatric and neurodegenerative disorders, as well as chronic pain. The complex structure of NMDARs permits several locations for therapeutic inhibition, making these receptors a potential target for multiple drugs which modulate them in different ways. NMDAR antagonists, which may be competitive, non-competitive, or uncompetitive, either block glutamate from binding the receptor or modulate the response to glutamate binding. Despite longstanding concerns about side effects of NMDAR antagonists, recent research suggests that, when appropriately used, these agents have favourable safety profiles. Furthermore, their fast-acting mechanism of action, resulting in rapid effects compared to other therapeutic agents, makes them a promising class of drugs that may yield effective therapeutics for multiple CNS disorders.

Contextual Constraints: Dynamic Evolution of Snake Venom Phospholipase A2.

Venom is a dynamic trait that has contributed to the success of numerous organismal lineages. Predominantly composed of proteins, these complex cocktails are deployed for predation and/or self-defence. Many non-toxic physiological proteins have been convergently and recurrently recruited by venomous animals into their toxin arsenal. Phospholipase A2 (PLA2) is one such protein and features in the venoms of many organisms across the animal kingdom, including snakes of the families Elapidae and Viperidae. Understanding the evolutionary history of this superfamily would therefore provide insight into the origin and diversification of venom toxins and the evolution of novelty more broadly. The literature is replete with studies that have identified diversifying selection as the sole influence on PLA2 evolution. However, these studies have largely neglected the structural/functional constraints on PLA2s, and the ecology and evolutionary histories of the diverse snake lineages that produce them. By considering these crucial factors and employing evolutionary analyses integrated with a schema for the classification of PLA2s, we uncovered lineage-specific differences in selection regimes. Thus, our work provides novel insights into the evolution of this major snake venom toxin superfamily and underscores the importance of considering the influence of evolutionary and ecological contexts on molecular evolution.

Role of Phospholipases A2 in Vascular Relaxation and Sympatholytic Effects of Five Australian Brown Snake, Pseudonaja spp., Venoms in Rat Isolated Tissues.

Human envenoming by Australian brown snakes (Pseudonaja spp.) may result in potentially life-threatening hypotension and subsequent cardiovascular collapse. There have been relatively few studies of the cardiovascular and sympathetic effects of Pseudonaja spp. venoms. In this study, we have examined the effects of venom from five brown snake species-P. affinis, aspidorhyncha, inframacula, nuchalis, and textilis-on cardiac inotropic and chronotropic responses, vascular tone, and sympathetic nerve-induced vascular contractions in rat isolated tissues. The role of phospholipases A2 (PLA2s) in venom-induced effects was assessed with the sPLA2 inhibitor varespladib. In rat isolated left and right atria, there were no physiologically relevant effects of Pseudonaja venoms (0.1-30 µg/ml) on left atrial force of contraction (inotropy) or right atrial rate (chronotropy). In contrast, in isolated small mesenteric arteries precontracted with a thromboxane mimetic, each of the five brown snake venoms (at 30 µg/ml) caused marked vasorelaxation (-60 to -90% of contractile tone). Pretreatment with varespladib (1 µM) significantly inhibited the vasorelaxation caused by P. aspidorhyncha, P. nuchalis, and P. textilis venoms. Electrically induced sympathetic nerve-mediated contractions of mesenteric arteries were significantly attenuated by only P. textilis, and P. affinis venoms (30 µg/ml) and these sympatholytic effects were inhibited by varespladib (1 µM). Based on their inhibition with the sPLA2 inhibitor varespladib, we conclude that PLA2 toxins in P. aspidorhyncha, P. nuchalis, and P. textilis venoms are involved in brown snake venom-induced vasorelaxation and the sympatholytic effects of P. affinis, and P. textilis venoms. Our study supports the promising potential role of varespladib as an initial (pre-referral) and/or adjunct (in combination with antivenom) therapeutic agent for brown snake envenoming.

Mutual enlightenment: A toolbox of concepts and methods for integrating evolutionary and clinical toxinology via snake venomics and the contextual stance.

Snakebite envenoming is a neglected tropical disease that may claim over 100,000 human lives annually worldwide. Snakebite occurs as the result of an interaction between a human and a snake that elicits either a defensive response from the snake or, more rarely, a feeding response as the result of mistaken identity. Snakebite envenoming is therefore a biological and, more specifically, an ecological problem. Snake venom itself is often described as a "cocktail", as it is a heterogenous mixture of molecules including the toxins (which are typically proteinaceous) responsible for the pathophysiological consequences of envenoming. The primary function of venom in snake ecology is pre-subjugation, with defensive deployment of the secretion typically considered a secondary function. The particular composition of any given venom cocktail is shaped by evolutionary forces that include phylogenetic constraints associated with the snake's lineage and adaptive responses to the snake's ecological context, including the taxa it preys upon and by which it is predated upon. In the present article, we describe how conceptual frameworks from ecology and evolutionary biology can enter into a mutually enlightening relationship with clinical toxinology by enabling the consideration of snakebite envenoming from an "ecological stance". We detail the insights that may emerge from such a perspective and highlight the ways in which the high-fidelity descriptive knowledge emerging from applications of -omics era technologies - "venomics" and "antivenomics" - can combine with evolutionary explanations to deliver a detailed understanding of this multifactorial health crisis.

How the Toxin got its Toxicity.

Venom systems are functional and ecological traits, typically used by one organism to subdue or deter another. A predominant subset of their constituent molecules-"toxins"-share this ecological function and are therefore molecules that mediate interactions between organisms. Such molecules have been referred to as "exochemicals." There has been debate within the field of toxinology concerning the evolutionary pathways leading to the "recruitment" of a gene product for a toxic role within venom. We review these discussions and the evidence interpreted in support of alternate pathways, along with many of the most popular models describing the origin of novel molecular functions in general. We note that such functions may arise with or without gene duplication occurring and are often the consequence of a gene product encountering a novel "environment," i.e., a range of novel partners for molecular interaction. After stressing the distinction between "activity" and "function," we describe in detail the results of a recent study which reconstructed the evolutionary history of a multigene family that has been recruited as a toxin and argue that these results indicate that a pluralistic approach to understanding the origin of novel functions is advantageous. This leads us to recommend that an expansive approach be taken to the definition of "neofunctionalization"-simply the origins of a novel molecular function by any process-and "recruitment"-the "weaponization" of a molecule via the acquisition of a toxic function in venom, by any process. Recruitment does not occur at the molecular level or even at the level of gene expression, but only when a confluence of factors results in the ecological deployment of a physiologically active molecule as a toxin. Subsequent to recruitment, the evolutionary regime of a gene family may shift into a more dynamic form of "birth-and-death." Thus, recruitment leads to a form of "downwards causation," in which a change at the ecological level at which whole organisms interact leads to a change in patterns of evolution at the genomic level.

Causes and Consequences of Snake Venom Variation.

Snake venoms are mixtures of toxins that vary extensively between and within snake species. This variability has serious consequences for the management of the world's 1.8 million annual snakebite victims. Advances in 'omic' technologies have empowered toxinologists to comprehensively characterize snake venom compositions, unravel the molecular mechanisms that underpin venom variation, and elucidate the ensuing functional consequences. In this review, we describe how such mechanistic processes have resulted in suites of toxin isoforms that cause diverse pathologies in human snakebite victims and we detail how variation in venom composition can result in treatment failure. Finally, we outline current therapeutic approaches designed to circumvent venom variation and deliver next-generation treatments for the world's most lethal neglected tropical disease.

Functional venomics of the Big-4 snakes of Pakistan.

In South Asia, the "Big-4" venomous snakes Naja naja, Bungarus caeruleus, Daboia russelii, and Echis carinatus are so-called because they are the most medically important snakes in the region. Antivenom is the only effective treatment option for snakebite envenoming but antivenom is not produced domestically in Pakistan making the country reliant on polyvalent products imported from India and Saudi Arabia. The present study investigated the toxin composition and activity of the venoms of Pakistani specimens by means of proteomic and physio/pharmacological experiments. To evaluate the composition of venoms, 1D/2D-PAGE of crude venoms and RP-HPLC followed by SDS-PAGE were performed. Enzymatic, hemolytic, coagulant and platelet aggregating activities of crude venoms were assayed and were concordant with expectations based on the abundance of protein species in each. Neutralization assays were performed using Bharat polyvalent antivenom (BPAV), a product raised against venoms from Big-4 specimens from southern India. BPAV exhibited cross-reactivity against the Pakistani venoms, however, neutralization of clinically relevant activities was variable and rarely complete. Cumulatively, the presented data not only highlight geographical variations present in the venoms of the Big-4 snakes of South Asia, but also demonstrate the neutralization potential of Indian polyvalent against the venom of Pakistani specimens. Given the partial neutralization observed, it is clear that whilst BPAV is a life-saving product in Pakistan, in future it is hoped that a region-specific product might be manufactured domestically, using venoms of local snakes in the immunising mixture.

Coagulotoxic effects by brown snake (Pseudonaja) and taipan (Oxyuranus) venoms, and the efficacy of a new antivenom.

Snakebite is a neglected tropical disease that disproportionately affects the poor. Antivenom is the only specific and effective treatment for snakebite, but its distribution is severely limited by several factors, including the prohibitive cost of some products. Papua New Guinea (PNG) is a snakebite hotspot but the high costs of Australian antivenoms (thousands of dollars per treatment) makes it unaffordable in PNG. A more economical taipan antivenom has recently been developed at the Instituto Clodomiro Picado (ICP) in Costa Rica for PNG and is currently undergoing clinical trials for the treatment of envenomations by coastal taipans (Oxyuranus scutellatus). In addition to potentially having the capacity to neutralise the effects of envenomations of non-PNG taipans, this antivenom may have the capacity to neutralise coagulotoxins in venom from closely related brown snakes (Pseudonaja spp.) also found in PNG. Consequently, we investigated the cross-reactivity of taipan antivenom across the venoms of all Oxyuranus and Pseudonaja species. In addition, to ascertain differences in venom biochemistry that influence variation in antivenom efficacy, we tested for relative cofactor dependence. We found that the new ICP taipan antivenom exhibited high selectivity for Oxyuranus venoms and only low to moderate cross-reactivity with any Pseudonaja venoms. Consistent with this genus level distinction in antivenom efficacy were fundamental differences in the venom biochemistry. Not only were the Pseudonaja venoms significantly more procoagulant, but they were also much less dependent upon the cofactors calcium and phospholipid. There was a strong correlation between antivenom efficacy, clotting time and cofactor dependence. This study sheds light on the structure-function relationships of the procoagulant toxins within these venoms and may have important clinical implications including for the design of next-generation antivenoms.

Proteomic and functional variation within black snake venoms (Elapidae: Pseudechis).

Pseudechis (black snakes) is an Australasian elapid snake genus that inhabits much of mainland Australia, with two representatives confined to Papua New Guinea. The present study is the first to analyse the venom of all 9 described Pseudechis species (plus one undescribed species) to investigate the evolution of venom composition and functional activity. Proteomic results demonstrated that the typical Pseudechis venom profile is dominated by phospholipase A2 toxins. Strong cytotoxicity was the dominant function for most species. P. porphyriacus, the most basal member of the genus, also exhibited the most divergent venom composition, being the only species with appreciable amounts of procoagulant toxins. The relatively high presence of factor Xa recovered in P. porphyriacus venom may be related to a predominantly amphibian diet. Results of this study provide important insights to guide future ecological and toxinological investigations.

Rattling the border wall: Pathophysiological implications of functional and proteomic venom variation between Mexican and US subspecies of the desert rattlesnake Crotalus scutulatus.

While some US populations of the Mohave rattlesnake (Crotalus scutulatus scutulatus) are infamous for being potently neurotoxic, the Mexican subspecies C. s. salvini (Huamantlan rattlesnake) has been largely unstudied beyond crude lethality testing upon mice. In this study we show that at least some populations of this snake are as potently neurotoxic as its northern cousin. Testing of the Mexican antivenom Antivipmyn showed a complete lack of neutralisation for the neurotoxic effects of C. s. salvini venom, while the neurotoxic effects of the US subspecies C. s. scutulatus were time-delayed but ultimately not eliminated. These results document unrecognised potent neurological effects of a Mexican snake and highlight the medical importance of this subspecies, a finding augmented by the ineffectiveness of the Antivipmyn antivenom. These results also influence our understanding of the venom evolution of Crotalus scutulatus, suggesting that neurotoxicity is the ancestral feature of this species, with the US populations which lack neurotoxicity being derived states.

Catch a tiger snake by its tail: Differential toxicity, co-factor dependence and antivenom efficacy in a procoagulant clade of Australian venomous snakes.

A paradigm of venom research is adaptive evolution of toxins as part of a predator-prey chemical arms race. This study examined differential co-factor dependence, variations relative to dietary preference, and the impact upon relative neutralisation by antivenom of the procoagulant toxins in the venoms of a clade of Australian snakes. All genera were characterised by venoms rich in factor Xa which act upon endogenous prothrombin. Examination of toxin sequences revealed an extraordinary level of conservation, which indicates that adaptive evolution is not a feature of this toxin type. Consistent with this, the venoms did not display differences on the plasma of different taxa. Examination of the prothrombin target revealed endogenous blood proteins are under extreme negative selection pressure for diversification, this in turn puts a strong negative selection pressure upon the toxins as sequence diversification could result in a drift away from the target. Thus this study reveals that adaptive evolution is not a consistent feature in toxin evolution in cases where the target is under negative selection pressure for diversification. Consistent with this high level of toxin conservation, the antivenom showed extremely high-levels of cross-reactivity. There was however a strong statistical correlation between relative degree of phospholipid-dependence and clotting time, with the least dependent venoms producing faster clotting times than the other venoms even in the presence of phospholipid. The results of this study are not only of interest to evolutionary and ecological disciplines, but also have implications for clinical toxinology.

Correlation between ontogenetic dietary shifts and venom variation in Australian brown snakes (Pseudonaja).

Venom is a key evolutionary trait, as evidenced by its widespread convergent evolution across the animal kingdom. In an escalating prey-predator arms race, venoms evolve rapidly to guarantee predatory or defensive success. Variation in venom composition is ubiquitous among snakes. Here, we tested variation in venom activity on substrates relevant to blood coagulation among Pseudonaja (brown snake) species, Australian elapids responsible for the majority of medically important human envenomations in Australia. A functional approach was employed to elucidate interspecific variation in venom activity in all nine currently recognised species of Pseudonaja. Fluorometric enzymatic activity assays were performed to test variation in whole venom procoagulant activity among species. Analyses confirmed the previously documented ontogenetic shift from non-coagulopathic venom in juveniles to coagulopathic venom as adults, except for the case of P. modesta, which retains non-coagulopathic venom as an adult. These shifts in venom activity correlate with documented ontogenetic shifts in diet among brown snakes from specialisation on reptilian prey as juveniles (and throughout the life cycle of P. modesta), to a more generalised diet in adults that includes mammals. The results of this study bring to light findings relevant to both clinical and evolutionary toxinology.